Cubic MnSb: epitaxial growth of a predicted room temperature half-metal 
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Epitaxial films including bulk-like cubic and wurtzite polymorphs of MnSb have been grown by 
molecular beam epitaxy on GaAs via careful control of the Sb4/Mn flux ratio. Nonzero-temperature 
density functional theory was used to predict ab initio the spin polarization as a function of reduced 
magnetization for the half-metals NiMnSb and cubic MnSb. In both cases, half-metallicity is lost 
at a threshold magnetization reduction, corresponding to a temperature T* < Tc- This threshold 
is far higher for cubic MnSb compared to NiMnSb and corresponds to T* > 350K, making epitaxial 
cubic MnSb a promising candidate for room temperature spin injection into semiconductors. 



Efficient electrical spin injection is crucial for the de- 
velopment of semiconductor (SC) spintronics, in which 
spin-polarized charge carriers flow from a ferromagnetic 
(FM) material into a non-magnetic SC structure. An 
ideal FM material would provide a fully spin polarized 
current [i] and should possess "engineering compatibil- 
ity" with its SC host, which includes controllable epitaxy 
and interfaces, and suitability for device processing. Cer- 
tain half-metallic ferromagnetic (HMF) materials have a 
density-of-states (DOS) spin polarisation at the Fermi 
level of Pdos = 100% due to the presence of a band 
gap straddling the Fermi energy Ep for one spin channel 
only [2*1 . Density functional theory (DFT) has predicted 
this property in several classes of materials [l| , including 
Heusler alloys and transition metal pnictides (TMPs) , at 
temperature T = OK. On this basis, half-metallicity was 
first envisaged in NiMnSb, which has a minority spin gap 
Eg PS 0.5 eV Q. However, there is presently no HMF can- 
didate with SC engineering compatibility which clearly 
maintains Pdos ~ 100% at room temperature. 

Experimental verification of Pdos is not straightfor- 
ward and after more than two decades' study, Pdos{T) 
for NiMnSb is not known. Positron annihilation exper- 
iments on bulk NiMnSb at T = 27K support Pdos ~ 
100% Q. Spin polarized photoemission spectroscopy 
(SPPES) can, in principle, measure Pdos more directly 
but it is strongly affected by the magnetic, chemical 
and structural conditions of the surface, which may dif- 
fer from those of the bulk. Threshold SPPES measured 
Pdos < 50% at T = 20K fH, while more recent works 
measured 40% [.6] or 50% at T = lOK falling to 40% 
at T = 300K 7i]. Inverse SPPES experiments implied 
similar Pdos values, again more dependent on surface 
condition than T [Sj. Note that spin polarizations in 
electron transport are not the same as Pdos, and have 
been inferred from tunnel junctions or spin light emit- 
ting diodes [l3] , or by Andreev reflection [ll| . For struc- 
tures incorporating NiMnSb, polarizations between 2% 
and 30% are typical, even at low temperatures. 

Given the difficulties of measuring Pdos unambigu- 
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ously, theoretical predictions play a key role in assessing 
the potential of HMF materials. The binary TMPs have 
attracted much recent attention. This class includes ma- 
terials such as MnAs and CrSb which are predicted to 
be wider-gap HMFs in cubic (sphalerite) [l3, and 
wurtzite [14,] (w-) structures. For example, we predict 
cubic c-MnSb to be HMF with Eg = 1.08 eV, while n- 
MnSb (its naturally occurring niccolite B81 structure) is 
a normal FM metal. Attempts have been made to sta- 
bilize cubic polymorphs by high-strain epitaxial growth 
on common SC substrates, e.g. c-CrAs and c-CrSb 
grown as ultra-thin films on GaAs(OOl). However, these 
films revert to their stable non-HMF structures within 
2-4 nm thickness [H [11] while X-ray diffraction (XRD) 
has shown the presence of strained orthorhombic o-CrAs 
rather than c-CrAs Tg']. DFT calculations predicted 
that biaxial strain would rarely stabilize such HMF poly- 
morphs [13] • Furthermore, a DFT study of CrAs ultra- 
thin films on GaAs(OOl) found that bulk-like o-CrAs was 
favorable for films with just 3 Cr atomic layers, and ar- 
gued that reported FM hysteresis is due to uncompen- 
sated interfacial spins of anti-FM o-CrAs 18|. One con- 
cludes that highly strained epitaxy is not a promising 
route to producing wide-gap HMF thin films of this class. 

While extrinsic effects can reduce spin polarization, 
such problems can be tackled by improving material qual- 
ity [TqI] and epitaxy [2^. It is crucial to address the in- 
trinsic limitations on Pdos- It has been realised that 
HMF band structures change as T increases, even well 
below the Curie temperature Tc [2l|-[2i]. The DOS in the 
minority spin gap becomes nonzero as the thermal fluc- 
tuations produce disorder in the local spin moment align- 
ment. In fact, beyond a threshold temperature T* < Tc, 
the value of Pdos collapses far more rapidly than the 
reduction of magnetisation M{T) [22! ]. 

We report here the successful stabilization of bulk- 
like c-MnSb and w-MnSb crystallites within fully re- 
laxed n-MnSb fllms grown by molecular beam epitaxy 
(MBE) on GaAs(lll). We confirm by DFT that these 
are HMF materials at T = OK. However, we also cal- 
culate Pdos as a function of magnetization M{T) using 
nonzero-temperature DFT, and show that Pdos in c- 
MnSb is much more robust against finite temperature 
effects than in NiMnSb. 



2 




FIG. 1. (color online) TEM data for polymorphic MnSb films 
on GaAs(fff). A typical bright-field image (a) and Fourier 
transforms from individual regions (b) are labelled with con- 
stituent phases, while a high resolution micrograph of the 
c-MnSb/n-MnSb boundary shows a sharp epitaxial interface. 

MnSb films were grown on ^10 mm wide GaAs(lll) 
substrates by MBE, without substrate rotation [25l[26|. 
The substrate temperature T^ub was fixed at 695 ± 5 K, 
and the growth rate was 6 nm min~^. For each growth 
run the Sb4/Mn fiux ratio J was measured directly by 
a retractable ionization gauge. Typically, samples were 
capped with Sb to prevent oxidation [23] • Specimens 
with varying total film thickness and J were prepared 
for transmission electron microscopy (TEM) using con- 
ventional techniques. They were examined using JEOL 
2000FX and Tecnai F20 microscopes operating at 200 kV, 
using bright field and dark field imaging as well as se- 
lected area diffraction patterns (SADP). 

In Fig. [l]we show TEM data from typical polymorphic 
films. Initially the epilayer grows as n-MnSb, with a high 
density of misfit dislocations at the n-MnSb/GaAs inter- 
face efficiently relaxing the 3.2% in-plane strain within a 
few nm thickness [2^,[2^ . However, beyond a certain film 
thickness (typically 100 nm) a sharp transition to poly- 
morph growth occurs. Careful analysis of the images and 
SADP reveals that two epitaxial polymorphs of MnSb ex- 
ist in the upper portion of the films, namely c-MnSb(lll) 
and w-MnSb(OOOl). Cahbrating SADP data to the GaAs 
substrate gives the following lattice parameters (in A, un- 
certainty ±0.5%): n-MnSb a=4.115, c=5.769; w-MnSb 
0=4.291, c=7.003; c-MnSb a=6.502. For n-MnSb, the 
lattice parameters conform with results in the literature, 
which gives us confidence in the w-MnSb and c-MnSb 



values, while for w-MnSb, the c/a ratio is almost exactly 
1.632 = y^8/3, corresponding to tetrahedral coordination 
of nearest neighbours. For the substrate and all three 
polymorphs, the half-crystal symmetries are 3m and in 
the absence of translational domains a single epitaxial 
orientation is predicted [2^. We see no TEM evidence 
for translational domains and associated double position- 
ing boundaries in the epilayer due to pre-existing GaAs 
surface steps. The polymorph growth is not influenced by 
post-growth cycling between 300 K and 700 K and the 
film structure is stable for many months, provided the 
samples are not heated outside ultra-high vacuum and 
are protected from oxidation Furthermore, XRD 

data (not shown) from as-grown films show the presence 
of w-MnSb and c-MnSb in addition to n-MnSb, and so we 
can rule out polymorphism as an artefact of TEM spec- 
imen preparation. The interface between n-MnSb and 
the polymorph epilayer is sharp, as shown by the high 
resolution image in Fig. [TJ 

On the basis of these observations, we attribute the nu- 
cleation of polymorphs to surface processes during MBE 
growth. In fact, the onset of polymorphism in the films 
is most strongly correlated with the fiux ratio J: careful 
control of both J and T^ub is required to produce c- and 
w-MnSb. For J < 5, the growth rate is Sb-limited [2^ 
while for J > 8 the film becomes faceted and non-planar. 
In the range 6.6 < J < 7.5, c-MnSb crystallites read- 
ily appear in the n-MnSb films. We have not observed 
polymorphism in films grown at T^ub < 675 K. Surface 
segregation of Sb leading to a double-layered surface re- 
construction with local tetrahedral coordination is a po- 
tential explanation for the transition, which requires a 
layer stacking order change from ABAC to AaBbCc or 
AaBbAa in the [0001] or [111] direction. The mecha- 
nism for stabilization of these polymorphs is clearly very 
different from direct epitaxial strain of the SC substrate 
previously investigated [H, [T5l - [l8j . XRD shows a pre- 
dominance of c-MnSb over w-MnSb in most films and we 
therefore focus on c-MnSb below. 

We calculated the T = electronic properties of 
the three polymorphs using two standard DFT packages 
[29I [soj . Full details will be given in a longer paper. 
Both c-MnSb and w-MnSb are HMF materials with Eg 
1.08 eV and 1.15 eV respectively and integer total spin 
moment AfiB- The calculated electronic structure of n- 
MnSb agrees with previous work [si'] and shows a weak 
Pdos of 18%. Trends in optimized DFT lattice parame- 
ters between the three polymorphs correspond well with 
the SADP data. The change of MnSb structure from 
n- to c- or w- represents true polymorphism since there 
is a change of nearest-neighbour bonding to tetrahedral. 
Conversely, the c- and w- phases are mutual polytypes, 
differing only in the stacking sequence of (111) or (0001) 
planes. As expected, they are very close in total energy, 
calculated respectively at 0.829 and 0.833 eV per formula 
unit above n-MnSb. 

We have extended the modelling of the c-MnSb and 
n-MnSb polymorphs to nonzero-temperature DFT us- 
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FIG. 2. (color online) Nonzero-temperature electronic struc- 
ture theory for MnSb and NiMnSb. The curves show the 
calculated Fermi level spin polarisation Pdos as a function of 
magnetization reduction for c-MnSb and NiMnSb, while the 
inset table summarises calculated Tc values from DLM-DFT. 
The experimental lattice parameter was assumed for c-MnSb, 
and the value used for NiMnSb matched that of Ref. [2211. 



ing the disordered local moment (DLM) theory [32!|, im- 
plemented using the Korringa, Kohn, Rostoker (KKR) 
multiple-scattering method |33l [33 | , to examine the ro- 
bustness of the half-metallicity. The theory includes 
the effects of thermally induced spin fluctuations on the 
underlying electronic structure, which must eventually 
destroy the long-range magnetic order and the overall 
spin polarization of the electronic states above the Curie 
temperature Tq- The DLM-DFT approach goes beyond 
the usual low temperature approximation of preset spin 
polarization and associates local spin polarization axes 
with each lattice site. These are allowed to vary slowly 
on the timescale of the electronic motions, affecting the 
nature and extent of the spin polarization of the elec- 
tronic structure and being self-consistently maintained. 
Thermally-induced disorder in the orientations of the lo- 
cal moments is treated using a mean- field scheme, based 
on the coherent potential approximation (35i . i36j . By 
taking appropriate ensemble averages over the local mo- 
ments' orientational configurations, we are able to de- 
termine the system's magnetic properties and find the 
temperature dependence of the spin polarized electronic 
structure and hence Pdos- This DLM-DFT theory also 
describes ah initio the onset and type of magnetic or- 
der [s^-fs^, l37l - [39| . The expression for the paramagnetic 

2 

spin susceptibility, x{^.^T) = 3fc^x„g(2)(q ; from which 
Tc values are estimated, is derived by considering the 
response of the DLM paramagnetic state to the applica- 
tion of an external, site-dependent, magnetic field. For- 
mally, 5*'^^ (q, T) is the direct correlation function for the 
local moments, of size /i, and is the analogue of an ex- 
change integral in the classical Heisenberg model of inter- 



acting spins. An important difference, however, is that 
5(^^(q, T) contains a fully itinerant description of elec- 
trons, whereas in the Heisenberg model the exchange in- 
tegral is defined in terms of interactions between localized 
spins. Details of the methods used to evaluate S^"^^ can 
be found in refs. [38l[39|. 

Results of our DLM-DFT calculations are summarised 
in Fig. [21 The inset tabulates the Curie temperatures 
for n-MnSb and NiMnSb: our calculated mean-field Tc 
values are larger than the experimental results but with 
the correct comparative magnitudes between the two ma- 
terials. Previous calculations 22| for NiMnSb have also 
somewhat overestimated Tc- Our calculated values for 
NiMnSb and c-MnSb are very similar. That Tc > 300 K 
for c-MnSb is promising but, as pointed out previously, 
does not guarantee a high Pdos at room temperature 
[13. The behavior of Pdos predicted by DLM-DFT as a 
function of magnetization reduction is plotted in Fig. ^ 
for c-MnSb and compared with NiMnSb: the horizontal 
axis is effectively a distorted temperature scale between 
and Tc- For NiMnSb, Pdos rapidly collapses even after 
around 5% reduction in magnetization. The same behav- 
ior is evident for c-MnSb but in this case Pdos remains 
at 100% to a much higher magnetization reduction of 
around 21%, essentially due to the much larger minority 
spin gap. Using DLM-DFT to estimate the magnetiza- 
tion dependence on temperature M{T) we can compare 
T* of NiMnSb and c-MnSb [13 • We find that while for 
NiMnSb T* > 100 K, in c-MnSb T* > 350 K, comfort- 
ably above room temperature. These values are lower- 
bound estimates since the theory typically overestimates 
the rate of decrease of M with rise in T; they are based 
on Tc « 730 K rather than the higher calculated values. 

In c-MnSb the nearest neighbours of each Mn atom are 
4 Sb atoms. The local tetrahedral environment enables 
the t2g states of Mn to hybridise easily with the Sb p- 
states so that a large bonding-antibonding splitting is 
set up where the lower-lying bonding states have more p- 
character whilst the higher antibonding states are more 
d-like. The gap in between is partly filled by the non- 
bonding narrow Eg states of the transition metal. Mn has 
a large exchange splitting which causes the bands of the 
majority and minority spin electrons to be distributed 
around different energies [i^. At T = the majority 
spin d-states of the Mn atoms are closer to the Sb p- 
states and the bonding-antibonding gap is smaller. For 
the minority spin electrons the Mn d-states are forced 
higher by the exchange splitting and the gap for this 
spin channel straddles Ep with the anti-bonding states 
unfilled. 

As the temperature is raised a local exchange splitting 
of the electronic structure persists even above Tc despite 
the orientations of the local spin polarization directions 
becoming disordered. A simple weighting of the majority 
and minority spin electronic structure to reflect the de- 
crease in overall spin polarisation would produce a linear 
dependence of Pdos on M{T). This is not seen in Fig. [21 
however, and to quantify T* we have to assess the effects 
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FIG. 3. (color online) Spin-polarized DOS for c-MnSb at 
different magnetizations M(T) calculated by DLM-DFT at 
the experimental lattice parameter. The three curves shown 
correspond to T = OK (blue), T = T* (red) and T ^ Tc 
(black). 



of the local moment disorder on the electronic structure 
more carefully. This is illustrated for c-MnSb in Fig. |3l 
The minority spin gap Eg — 1.08 eV appears in the fully 
magnetized case corresponding to T = OK (blue curve). 
The separation between the antibonding t2g states of the 
majority and minority spin channels is large so that as the 
overall magnetization is reduced, "impurity" states asso- 
ciated with minority spin electrons appear in the major- 
ity spin channel above Ep whilst majority spin impurity 
states form in the minority spin gap close to the posi- 
tion of the antibonding t2g majority spin band. As M 
is reduced, these encroach up to the Fermi level and at 
T* remove the half-metallicity abruptly. This is shown 
in the red curve of Fig. |31 the impurity states have just 
reached Ep at a magnetization reduction of 21%, initi- 
ating the collapse of Pdos shown in Fig. [5] Finally, for 
T = Tc all spin polarization is lost (black curve). For 



NiMnSb, the same effects occur but T* corresponds to 
a magnetization reduction of just 5% due mainly to the 
smaller value of Eg . 

The mechanism of spin polarization reduction in 
hybridization-gap materials has been confirmed by the 
present fully ah initio calculations and elsewhere |2l| - [24| . 
It is clear that, as well as Tc, Eg is a key parameter 
in determining the high-temperature spin polarization. 
Gaining full control of the MBE growth of TMP poly- 
morphs with large gap will be an essential step towards 
exploiting these new materials. Fully in situ investiga- 
tions of the transition to polymorphic growth should shed 
more light on the microscopic surface processes leading 
to the growth transition. The combination of engineer- 
ing compatibility, large minority spin gap and robust spin 
polarisation make c-MnSb highly attractive for spintron- 
ics applications. Epitaxial interfacing to a semiconduc- 
tor remains a crucial task but the prognosis here is also 
good. MoUet and Jenkins [4l[ have calculated the (001) 
and (111) surface electronic structure of c-MnSb using 
DFT and found that half-metallicity is preserved for the 
favoured Sb-terminated surfaces. Since this structure is 
shared with common semiconductors such as GaAs, it is 
likely that half-metallic interfaces can readily be formed. 
In contrast, the preservation of half metallicity at SC in- 
terfaces is much less common with ternary Heusler alloys 
[20| . Interestingly, Jenkins's recent prediction of tilted 
Dirac cones (similar to those found in graphene and at 
topological insulator surfaces) in the surface states of 
NiMnSb [i^ is also likely to translate to c-MnSb, since 
the relevant symmetry of the structures is identical. 
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